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We study the non-Markovian entanglement dynamics of two qubits in a common squeezed bath. 
We see remarkable difference between the non-Markovian entanglement dynamics with its Markovian 
counterpart. We show that a non-Markovian decoherence free state is also decoherence free in the 
Markovian regime, but all the Markovian decoherence free states are not necessarily decoherence free 
in the non-Markovian domain. We extend our calculation from squeezed vacuum bath to squeezed 
thermal bath, where we see the effect of finite bath temperatures on the entanglement dynamics. 
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I. INTRODUCTION 

Entanglement is a remarkable feature of quantum me- 
chanics, and its investigation is both of practical and 
theoretical significance. It is also viewed as a basic 
resource for quantum information processing [l|, like 
entanglement-based quantum cryptography [2| , quantum 
telcportation Q , dense coding Q and cluster-state quan- 
tum computation Q. Thus the real world success of 
these quantum information processing schemes relies on 
the longevity of entanglement in multiparticle quantum 
states. Entanglement is also related to the basic issue of 
understanding the nature of nonlocality in quantum me- 
chanics p, |7|- However, a quantum system used in quan- 
tum information processing inevitably interacts with the 
surrounding environment, which induces the quantum 
world into the classical world [g, 3 ■ The presence of de- 
coherence in communication channels and computing de- 
vices degrades the entanglement when the particles prop- 
agate or computation evolves. The coupling of the quan- 
tum system with its surroundings and the consequent de- 
cay of entanglement motivate important questions such 
as to understand its sources, and possibly to find ways 
to circumvent it through different types of controlled en- 
vironments. 

The dynamics of open quantum systems, however, may 
be rather involved, mostly due to the complex structure 
of the environment interacting with the quantum system. 
Generally, the nonunitary evolution of the reduced den- 
sity matrix of the system is obtained after taking partial 
trace over bath variables. In this process, some approx- 
imations are often made in the derivation of a master 
equation for the systems reduced density matrix. The 
most important approximations [8|, [l0( are the weak cou- 
pling or Born approximation assuming that the coupling 
between the system and the reservoir is small enough 
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to justify a perturbative approach, and the Markov ap- 
proximation assuming that the correlation time of the 
reservoir is very short compared to the typical system 
response time so that the reservoir correlation function 
is assumed to be ^-correlated in time. Although, the use 
of Markovian approximation is justified in a large variety 
of quantum optical experiments where entanglement has 
been produced, one should notice that non-Markovian ef- 
fects are crucial, e.g., for high-speed quantum communi- 
cation where the characteristic time of the relevant sys- 
tem become comparable with the reservoir correlation 
time, or if the environment is structured with a partic- 
ular spectral density, e.g. for quantum systems (chan- 
nels) embedded in solid-state devices, where memory ef- 
fects are typically non-negligible. In these cases, the dy- 
namics can be substantially different from the Markovian 
one. Due to their fundamental importance in quantum 
information processing and quantum computation, non- 
Markovian quantum dissipative svstems have attracted 
much attention in recent years |8|, Illl - l22f , one of the main 
purposes of which in the long run is to engineer different 
types of (artificial) reservoirs, and couple them to the sys- 
tem in a controlled way. The non-Markovian features of 
system-reservoir interactions have made great progress, 
but the theory is far from completion, in particular how 
different kinds of non-Markovian environments influence 
the systems, and the difference between Markovian and 
non-Markovian system evolutions are still a subject that 
demands further investigations. 

If the environment would act on the various parties 
the same way it acts on single systems, one would ex- 
pect that a measure of entanglement, say the concur- 
rence, would also decay exponentially in time. However, 
this is not always the case. Recently Yu and Eberly 
|23( showed that under certain conditions, the dynam- 
ics could be completely different and the quantum en- 
tanglement of a bipartite qubit system may vanish in a 
finite time. They called this effect "entanglement sudden 
death (ESD)". This phenomenon of entanglement sudden 
death has been extensively studied [2J] in the context of 



Markovian master equation. Entaglement dynamics for 
system of two qubits interacting with their local indepen- 
dent reservoir is completely different from that when they 
interact with the same common reservoir. We see that 
entanglement may not be destroyed by the interaction 
with the environment, and sometimes it persists at what- 
ever be the temperature of the bath. A common reser- 
voir indirectly couples the qubits, and there have been 
some suggestions |25j for creating as well as enhancing 
|26j entanglement between two or more parties by their 
collective interactions with a common environment. In 
the same vein, non-Markovian entanglement dynamics is 
fundamentally different from the Markovian one. With 
respect to the study of non-Markovian dynamics of two 
independent qubits, each locally interacting with its own 
reservoir, it was shown [18|, [l9[ that although no interac- 
tion is present or mediated between the qubits, there is 
a revival of their entanglement following the ESD. The 
backaction of the non-Markovian reservoir [18|, [l9| is re- 
sponsible for revivals of entanglement after sudden death. 
The dynamics of entanglement in two independent non- 
Markovian channels is shown [2(| to be oscillating at high 
temperatures, whereas in the Markovian channel, entan- 
glement was shown to decay exponentially. Appearance 
of sudden death and sudden birth of entanglement was 
also discussed [2l| in common structured reservoirs. 

Different schemes have been derived to remove the ef- 
fects produced by the environment, for example, quan- 
tum error correction, decoherence-free subspace (DFS), 
dynamical decoupling and quantum Zcno effect. In the 
presence of the environment, the DFS is a set of all states 
which is not affected at all by the interaction with the 
bath. There were some proposals related to the use of 
DFS as the memory space for storing the quantum in- 
formation [23| . Recently, decoherence free entanglement 
was studied for two two-level systems interacting with 
a common squeezed vacuum bath [28|. The Markovian 
entanglement dynamics of two two-level atoms that in- 
teract with a common squeezed vacuum reservoir was ex- 
tensively studied [22l ■ The phenomenon of sudden death 
and revival of entanglement was investigated for the ini- 
tial states that are very close to (as well as far from) the 
Markovian DFS. It was claimed that for states belonging 
initially to the DFS plane, the phenomenon of entangle- 
ment sudden death never occurs. However, if the ini- 
tial state is away from the DFS plane, the sudden death 
shows up, followed by sudden revival of entanglement. 

Several proposals to physically realize the squeezed 
reservoir were put forward in the literature, the sim- 
plest of which consists in considering a two-level atom 
immersed in a squeezed multimode radiation field |32| — 
|34| . Parkins et al. [35[ showed how a two-level system 
can be coupled to an almost ideal squeezed vacuum by 
assuming an atom strongly interacting with a cavity field 
which is illuminated by finite-bandwidth squeezed light. 
It was also shown |37| how a squeezed environment can 
be obtained by means of a suitable feedback of the out- 
put signal corresponding to a quantum-nondcmolition 



(QND) measurement of suitable quadrature operators. 
In Ref. [36] , the authors showed how to mimic the in- 
teraction of a two-level system with a squeezed reservoir 
by using a four-level atom interacting with circularly po- 
larized laser fields. Assuming a strong decay of the two 
most excited levels, it was shown that the dynamics of 
the two ground atomic states is effectively similar to that 
of a two-level system interacting with a squeezed reser- 
voir. A recent theoretical study [38( was made to generate 
squeezed reservoir for a two-level system by engineering 
the Hamiltonian of a A-type three-level atom interacting 
with a single cavity mode and laser fields with suitable 
intensity and detuning. Another experimental proposal 
[39t | for two two-state atoms in a common squeezed reser- 
voir was made employing quantum-reservoir engineering 
to controllably entangle the internal states of two atoms 
trapped in a high-finesse optical cavity. Using laser and 
cavity fields to drive two separate Raman transitions be- 
tween stable atomic ground states, a system, correspond- 
ing to a pair of two-state atoms coupled collectively to a 
squeezed reservoir, could be realized. 

In this paper, we will analyze the entanglement dynam- 
ics for two qubits interacting with a common squeezed 
reservoir in the non-Markovian regime. We see multiple 
cycles of entanglement sudden death and revival in the 
non-Markovian case, showing striking difference between 
the Markovian and non-Markovian entanglement dynam- 
ics. We extend our result for the finite temperature case 
where we see the non-Markovian entanglement oscilla- 
tions gradually decreases as one increases the tempera- 
ture. Finite temperature of the bath accelerates the phe- 
nomenon of ESD in general. We show that the Marko- 
vian decoherence free states remains invariant under fi- 
nite temperature of the bath, and all states in the Marko- 
vian DFS plane is not necessarily decoherence free in the 
non-Markovian regime. Interestingly, the singlet state 
(which satisfy a more general DFS condition) is found to 
be decoherence free both in the Markov and non-Markov 
regime, and is also found to be robust against finite bath 
temperature. In Sec. UTJ we describe our model and 
present the non-Markovian master equation. In Sec. IIIII 
we discuss the difference of DFS and entanglement dy- 
namics between the Markovian and non-Markovian cases. 
In Sec. IIVI we present the numerical results and discuss 
our main observations, ending it with some concluding 
remarks in Sec. [Vl 



II. MODEL AND QUANTUM MASTER 
EQUATIONS 

We consider a pair of two-level atoms (two qubits) 
coupled to a common non-Markovian thermal squeezed 
reservoir. The microscopic Hamiltonian of the system 
plus reservoir is given by 



H = ^ {a 1 +( j 1 _+a 2 + a 2 _)+J2^kb k t b k + H I , (1) 



where the interaction Hamiltonian Hj has the form 

H I = Y / 9kS+b k +g* k S-b k l 



(2) 



Here S + = a\_ + a\ and S- = a\_ + <j 2 _ are collective rais- 
ing and lowering operators for the two-qubit system with 
cr+ = |li)(0,-|, &L = |0i)(li|, where |lj) and |0j) are up 
and down states of the ith qubit, respectively. Let us now 
proceed for the master equation for this two-qubit sys- 
tem interacting with a common squeezed thermal bath 
according to the Hamiltonian (J]) . We assume the factor- 
ized initial system-reservoir state with the initial state 
of the reservoir as a squeezed thermal equilibrium state 
given by 
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(5) 



where f3 = 1/KT with K being the Boltzman constant 
and T being the temperature, and we have introduced the 
unitary squeeze [8| operator U{r k ,9 k ) with £/. = r k e l6k . 

The non-Markovian master equation in the interaction 
picture for the two-qubit reduced density matrix p{t) in 
the Born approximation can be calculated and written 
as: 



dp 

dt 



A(i) {S+pS- - pSS+} + A*(t){S+pS- - SS+p} 

+ p{t) {S- P S+ - S+S-p] + p*(t) {S- P S+ - pS+S-} 

+ a(t) {2S+pS+ — S+S+p — pS+S+} 

+ a*(t) {2S_ P S- - SS-p - pS^S-}. (6) 



The time dependent coefficients appearing in the master 
equation are respectively given by 

A(t) = / dti / duj J(u>)N(u;)e i ( Uo - u K t ~ tl \ (7) 



H(t)= / dh / dio J(w)[l + ^(w)]e 4(wo ^ )(t -* l) (8) 



a(t) 



where 



dt-y 



dto J{u)M{w)e i{ - u °- u){t+tl \ 



(9) 



N(u)=n(u) [cosh 2 r + sinh 2 r] + sinh 2 r, (10) 
M(u) = - coshr sinhr e lB [2n(w) + 1] (11) 

with the frequency independent resonant squeeze pa- 
rameter r and the resonant phase 9, and we have 



written n(uj) = 1/ [exp(/3w) — 1] for the Planck dis- 
tribution. The master equation ([5]) is valid for ar- 
bitrary temperature (provided that the Born approxi- 
mation still holds) and the squeezed vacuum reservoir 
case is just the zero-temperature limit of it. In the 
zero-temperature limit, n(ui) = 0, so, N(w) = N = 
sinh 2 (r), and M(uj) = — cosh(r) sinh(r) e l °— —M e l$ , 
where M = sinh(r) cosh(r) = ^N(N + 1). The non- 
Markovian character is contained in the time-dependent 
coefficients, which contain the information about the 
system-reservoir correlations. In the previous equations, 
J(uj) =J2k\9k\ 2 3( UJ — uj k ) is the spectral density charac- 
terizing the bath, where the index k labels the different 
field mode of the reservoir with frequency u> k . We may 
consider any form of the reservoir spectral density. But 
for simplicity, here we consider an Ohmic squeezed bath 
[ill fl2j with the spectral density given by 



J(w) = Tuj exp I 



2 K 



(12) 



where uj is the frequency of the bath and w c is the high- 
frequency cutoff and T is a dimcnsionless constant char- 
acterizing the interaction strength to the environment. 
For finite temperature Markovian case, n(ui) = n(w ). 
Since we consider two qubits in a common bath, they 
will have to be quite near so that the interatomic sepa- 
ration is much smaller than a typical wavelength of the 
bath, i.e., the length scale of the resonant wavelength 
Ao = hc/ojQ in the model [40| , where c is the wave speed of 
the bath. Thus we do not consider here the effect of qubit 
size, also the spatial dependence on qubit-environment 
coupling strength [40, HJJ {g k is assumed to be position- 
independent). We also assume that there is no direct 
interaction between the qubits except the indirect cou- 
pling through the common environment. 



III. DIFFERENCE BETWEEN MARKOVIAN 
AND NON-MARKOVIAN DYNAMICS 

Our next aim is to show the difference between the 
Markovian and non-Markovian entanglement dynamics 
for this system-reservoir model. The Markovian mas- 
ter equation, in the interaction picture, for two two-level 
systems interacting with a broadband squeezed vacuum 
bath (at zero temperature) is well studied |28l - l30l ] and is 
given by 



dp 

In 



~f{N + 1) (2S^pS+ - S+S^p - pS+S_) 
7-/V (2S + pS- — S-S + p — pS-S + ) 
1 Me ie (2S +P S+ - S + S+p - pS+S+) 



- -jMe" w (2S-_p5_ - 5_5_/9 - pSS-.) (13) 

where 7 = 7rru;o- Using this Markovian master equa- 
tion (|13[) . the entanglement dynamics (the phenomenon 
of sudden death and revival of entanglement) of a pair of 



two- level atoms has been extensively studied [29(. The 
DFS for this model (Markov approximation, broadband 
squeezed vacuum) is found in Ref. [28|, and we call it 
Markovian DFS. 

The main result of the present paper will rotate around 
the discussion of DFS and entanglement dynamics ac- 
cording to the Markovian master equation (|13[) and its 
non-Markovian counterpart ([5]), showing striking differ- 
ence between them. This Markovian master equation 
(|13[) can be written in an explicit Lindblad form using 
only one Lindblad operator [281 ] : 



dp 

at 
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L = VN + 1 S- 



N cxp{i6>} £4 



(14) 



(15) 



In this case, the DFS [42J is composed of all eigenstates 
of L with zero eigenvalues. The two orthogonal vectors 
in the DFS plane for this Markovian evolution are [28[ 
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One can also define the states ^3) and \cf>4) orthogonal 
to {|0i),|(/>2)} plane: 



|0 3 ) = - 7 =(|O1> + |1O)), 



(18) 
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It is important to mention here that the state |</>i), 
Eq. (fTS)) . is a decoherence free state only for the Mako- 
vian evolution (TLSf , but for a general non-Markovian dy- 
namics \4>\) is not decoherence free. A more general dis- 
cussions on DFS condition was discussed in Refs.j42l.l43l]. 
It is important to note that decoherence is the result of 
the entanglement between system and bath caused by the 
interaction term Hi , Eq. @ , of the Hamiltonian (fTJ . In 
other words, if Hi = then system and bath are decou- 
pled and evolve independently and unitarily under their 
respective Hamiltonians Hg and Hg. Clearly, then, a 
sufficient condition for decoherence free dynamics is that 
Hj = 0. However, since one cannot simply switch off 
the system-bath interaction, in order to satisfy this con- 
dition, it is necessary to look for special subspace (say, 
H) of the full system Hilbcrt space such that the system 
evolves in a completely unitary fashion on that subspace 
H. As shown first by Zanardi and Rasetti [43], such a 
subspace is found by assuming that there exists a set of 
degenerate eigenvectors of the system coupling operators 
S± in the system-reservoir interaction Hamiltonian. In 
our case, focusing on the form of interaction Hamiltonian 
Hj given by Eq. ©, it is clear that the DFS is made up 
of those states \ip) satisfying [42| 



The singlet state, Eq. ([17)) . is one which satisfies ([20)) 
with vanishing total angular momentum. Hence the sin- 
glet state \ip) = 1^2) = (|01) — |10))/-\/2 is a decoherence 
free state for this type of Hamiltonian. This is why the 
singlet state \4>2) is decoherence free both for Markovian 
([131) and non-Markovian ([5]) dynamics at any temper- 
ature. This result is also confirmed by the numerical 
calculations shown in Sec. IIV1 

We note here that the qualitative characteristics of 
the non-Markovian entanglement dynamics with oscilla- 
tions and sudden deaths and revivals (shown in Sec. IIV[) 
may not be specific to the squeezed reservoir, but the 
quantitative feature of this dynamics will depend on var- 
ious properties of the reservoir, such as spectral density, 
squeezing, and temperature. Nevertheless, it is impor- 
tant to mention here that the Markovian master equa- 
tion for two two-level system interacting with a common 
(unsqueezed) heat bath is 

% = \l (N + 1) (2S^ P S+ - S + S-p - pS+S^) 

+ -7JV (2S+pS- - SS+p - pSS+) . (21) 

In this case, the DFS is composed of common eigenstates 
of the Lindblad operators L± ~ S± with zero eigenvalues. 
The singlet state |^ 2 ) = (|01) — |10))/v / 2 satisfies this con- 
dition. So, the singlet state Ifo) is the only state which 
is decoherence free for an unsqueezed common reservoir. 
Interestingly, we see that for an unsqueezed reservoir, 
DFS calculated from the Lindblad operators L± = S± 
and that obtained from the system-bath interaction (Za- 
nardi and Rasetti criterion [42, |43|) Hamiltonian Hi, arc 
the same. Hence there will be no difference between the 
Markovian and non-Markovian entanglement dynamics 
for this DFS state. On the other hand, for the squeezed 
reservoir, one can have infinitely many Markovian DFS 
states |0i), Eq. (flo) , just by varying continuously the 
squeeze parameters 9 and r. In this case, squeezing 
plays an important role in showing difference between the 
Markovian and non-Markovian entanglement dynamics 
for these Markovian DFS states. This is quantitatively 
shown in Fig. [3] in the next section. Now, in order to 
study the sudden death and revival of entanglement of 
two qubits in this common squeezed bath (both in the 
Markovian and non-Markovian regimes), we consider as 
initial states of the form [H [H 



|vI/ 1 )= e | ( /, 1 ) + V /l- e 2| ( /, 4 ), 



|*2>=e|02> + Vl-e 2 |03), 



(22) 



(23) 



s±W = o. 



(20) 



where e is a variable amplitude of one of the states be- 
longing to the Markovian DFS plane {|0i), |02)}- We 
would like to study the effect of varying e (0 < e < 1) on 
the sudden death and revival of entanglement for these 
initial states. We calculate the time-evolved two-qubit 
density matrix for the initial states |\E'i) and l^) us- 
ing the Markovian ([13)) as well as the non-Markovian ^ 
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FIG. 1. (Color online) Non-Markovian (solid lines) and 
Markovian (dashed lines) time evolution of the concurrence 
for |\Pi) as an initial state with varying e: (a) e = 0, (b) 
e = 0.5, (c) e = 0.9, and (d) e = 1. The values of the other 
parameters are r = 0.31, = 0, 7 = 1, and ui c — uiq = 1. 
The number of ESD and revivals in the non-Markovian case 
increases. The state |</>i) (|^i) with e = 1) no longer remains 
decoherence free for the non-Markovian evolution although it 
is decoherence free in the Markovian case. 



master equations. The various components of the time- 
dependent density matrix depend on the initial states as 
well as on the squeezing parameters. 



IV. NUMERICAL RESULTS 

We calculate numerically the time evolution of the den- 
sity matrix according to Eqs. (|L3|) and® and we choose 
the Wootters entanglement measure |44j , the concurrence 
C(t), defined for the time-evolved two-qubit density ma- 
trix p(t) as 



C(t) = max(0, VAi - V^2 - v% - y/X*), (24) 

where Ai, A2, A3, A4 are the eigenvalues of the matrix 
p c = p{t){<jy <g> <jy)p*(t)(a y ® a y) in descending order. 
The entanglement dynamics (in both the Markovian and 
non-Markovian regimes) is shown in Fig.Q]for two qubits 
initially in the state |^i). The state |^i) is a superposi- 
tion of two states \<j>{) (belonging to the Markovian DFS) 
and its orthogonal \<j>i). We vary e between and 1 for 
fixed values of the parameters r = 0.31 and 9 = as 
in RefJ29[ . Let us recapitulate the observation made in 
Ref. [29[ about the Markovian dynamics of entanglement 
for the state |^i) in the parameter interval < e < 0.5 
[see also the curve in dashed line in Fig. Ufa)], where it 
was shown that the initial entanglement decays to zero 
in a finite time td- Then after a finite period of time 
during which the concurrence stays null, it revives at a 
later time t r , then reaching asymptotically its steady- 
state value. For this Markovian case, it was also observed 
that this death and revival cycle happens only once for 
the initial state I'J'i). For e = 0.5, the entanglement 
dies and revives simultaneously and eventually goes to 



its steady-state value [see also the curve in dashed line 
in Fig. [lfb)]. For 0.5 < e < 1, no entanglement sudden 
death was found [29| in the Markovian dynamics [see also 
the curve in dashed line in Fig.QJc)]. Finally, when e = 1, 
l^i) = 1 0i ) is a decoherence- free state so that the con- 
currence remains constant with time [see also the curve 
in dashed line in Fig. (Hd)]. This was the picture for the 
Markovian master equation (|13|) . Now we compare this 
Markovian case with the non-Markovian entanglement 
dynamics according to Eq.(j6|) for the initial state |^i). 

In the non-Markovian case for e = 0, we see from the 
curve in solid line in FigOJa) that the initial entangle- 
ment decays to zero in a finite time showing suddent 
death of entanglement and then it revives again, following 
four successive death and revival, and finally the steady 
state value of the concurrence is reached at large times. 
Whereas in the case of Markovian dynamics, we see only 
one death and revival. It is also important to note that 
the entanglement sudden death occurs much faster in 
the non-Markovian case compared to the Markovian one. 
The time gap between adjacent death and revival is small 
(that is the rate at which the death and revival occurs are 
very fast) compared to its Markovian counterpart. The 
striking difference shown in FigQJa) between the Markov 
and non-Markov entanglement dynamics is that in the 
non-Markovian case, the entanglement is nonzero (show- 
ing three revival cycle) in a time window when the con- 
currence remains null in the case of Markovian dynamics. 
Similar multiple death and revival cycle is observed for 
the state with e = 0.5 in the non-Markovian case showing 
a clear departure from the Markovian dynamics. When 
0.5 < e < 1, that is when we get closer to the Markovian 
DFS, it was shown [29J that the whole phenomenon of 
sudden death and revival disappears for the initial state 
l^i). Contrary to that, we see [for the initial state |^i) 
with e = 0.9 in Figfljc)] clear sudden death and revival 
for this range of e as well in the non-Markovian case. 

Fig. [2] shows the dynamical behavior of the entangle- 
ment in terms of the concurrence for the initial state | "J/2) - 
The state l^) is a superposition of two states 1^2) (be- 
longing to the Markovian DFS) and its orthogonal \</>3). 
We again vary e between and 1 for fixed values of the 
parameters r = 0.31 and 9 = 0. For the initial state 
1^2), we also see multiple cycles of death and revival 
of entanglement in the non-Markovian regime (in the 
parameter interval < e < l/v2) showing completely 
different behavior from its Markovian counterpart [29j |. 
When l/-\/2 < e < 1, no sudden death was observed for 
the state | \I/ 2) in the Markov regime whereas we see clear 
sudden death in this case also for the non-Markovian case 
[see Fig. Eld)]. 

Finally we emphasize again that the state \4>i) (the 
state |^i) with e = 1) is not a decoherence free state 
for non-Markovian master equation ^ although it be- 
longs to the DFS for Markovian master equation (fT~3]) 
[see Fig. Did)]. The Markovian and non-Markovian en- 
tanglement dynamics for the Markovian DFS states \4>i), 
Eq. (|16[) . with varying squeeze parameters r and 9 is 
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FIG. 2. (Color online) Non-Markovian (solid lines) and 
Markovian (dashed lines) time evolution of the concurrence 
for | $2) as an initial state with varying e: (a) e = 0.1, (b) 
e = 0.4, (c) e = 0.54, and (d) e = 0.707. The values of the 
other parameters are r = 0.31, 9 — 0, 7 = 1, and uj c = luq = 1. 
The number of ESD and revivals in the non-Markovian case 
increases. The state \<j>2) (1*2) with e = 1) remains decoher- 
ence free for both the Markovian and non-Markovian evolu- 
tion. 
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FIG. 3. (Color online) Non-Markovian entanglement dy- 
namics for various Markovian DFS states \<pi) with varying 
squeeze parameters r and 9: (a) r — 0.05 and r — 0.09 with 
= 0, (b) 9 = 7r/6 and = tt with r = 0.3. The values of the 
other parameters are 7 = 1, and w c = Wo = 1- 



shown in Fig. [3] Wc consider in Fig. [3] two separate cases 
(a) r = 0.05, r = 0.09 with 9 = and (b) (9 = tt/6, 
9 = tt with r = 0.3. We find that the initial entangle- 
ment and final asymptotic entanglement at large times 
(Markov and non-Markov) for the initial Markovian DFS 
states, Eq. (|T^|) . depend on r not on 9. One can ob- 
serve from Fig. [3^b) that when the squeeze parameter r 
is fixed and the phase 9 is varied, the Markovian concur- 
rence curves do not vary and they remain the same and 
overlap. This is because when r is varied, the values of 
N and M in Eq. (fT6]) change considerably and so do the 
resultant Markovian DFS states, while different values of 
9 change only the relative phase between |11) and |00) in 
the resultant Markovian DFS states of Eq. p6[) . Never- 
theless, one can see that these Markovian DFS states are 
not decoherence-free in the non-Markovian regime and 
the characteristics of this non-Markovian oscillations of 
concurrence depend on the squeeze parameters r and 9. 



FIG. 4. (Color online) Non-Markovian and Markovian time 
evolution of the concurrence at finite temperatures for |\Pi) as 
an initial state with (a) e — and (b) e — 1, and for 1*2 ) as an 
initial state with (c) e = 0.1 and (d) e = 0.54. The values of 
the other parametrs are r = 0.31, 9 — 0, 7 = 1, and u c = loq = 
1. The insets show the entanglement dynamics in the short- 
time region. The non-Markovian entanglement oscillations 
gradually disappear as one increases the temperature. The 
state |(jf>2) remains decoherence free both in the Markov and 
non-Markov regime at any temperature of the reservoir. But 
the state |<^i) (|*i) with e = 1) that is decoherence free at 
any reservoir temperature in the Markovian case no longer 
remains decoherence free for the non-Markovian evolution. 



On the other hand, the state |</>2) (the state 1*2) with 
e = 1) is decoherence free both for Markovian and non- 
Markovian master equations. This was verified by ob- 
taining a straight line C(t) = 1, showing that the entan- 
glement remains constant for this state at zero tempera- 
ture of the reservoir. 

Next, we go to the finite-temperature case for which 
N(u) and M(u) are given by Eqs.flTDJl and (JTTJ- To show 
the effect of temperature on the Markovian evolution let 
us focus on the entanglement dynamics of two specific 
initial states (|'Fi) with e = 0, and 1*2) with e = 0.1). 
From the zero-temperature Markovian dynamics of these 
two states [see Fig.QJa) and Fig.^a)] we see that at long 
time the entanglement finally saturates to a finite value 
after which no death of entanglement occurs, whereas at 
finite temperatures (say, KT > 2uj$) Markovian dynam- 
ics of these two states show a complete death of entan- 
glement at long-time limit [see Figs. HJa), Hfc) and their 
insets] . From Fig. [2ja) we see for the zero-temperature 
case, Markovian entanglement sudden death occurs at 
relatively slower rate whereas in the finite-temperature 
case [Fig. EJc)] ESD rate is much faster. The differ- 
ence between the zero-temperature entanglement dynam- 
ics and the finite-temperature entanglement dynamics is 
more transparent in the non-Markovian case. The non- 
Markovian entanglement dynamics at zero temperature 
shows several death and revival cycles in short time limit 
while there is only one or two death-revival cycles fol- 
lowing a decay of entanglement for the finite tempera- 
ture (say, KT > 2ujq) non-Markovian case. In the non- 



Markovian evolution (Fig. [4j, as one gradually increases 
the temperature, the number of death and revival cycles 
of the concurrence reduces and the effect of temperature 
on concurrence in general is seen as decaying in nature 
(concurrence asymptotically reaches to its zero value at 
large times) with an exception to the DFS. That is in the 
long time limit, the non-Markovian entanglement satu- 
rates to a certain nonzero value at low temperature, while 
at high temperatures {KT > 2wo) the concurrence finally 
decays to zero in the long-time limit. Note also that the 
concurrence disappear more quickly at finite tempera- 
tures for both Markovian and non-Markovian cases. It 
is seen that the death time for the Markovian dynamics 
is always less compared to the non-Markovian dynam- 
ics at a given finite temperature. We have seen that the 
state |</>2) (the state |\&2) with e = 1) remains decoher- 
ence free both in the Markov and non-Markov regime at 
any finite temperature of the reservoir. This was verified 
by obtaining a straight line C(i) = 1, showing that the 
entanglement remains constant for this state \cf>2) at any 
arbitrary temperature of the reservoir. It is quite striking 
that the state \<f>i) shows entanglement sudden death at 
finite temperature in the non-Markovian dynamics, al- 
thouth it remains a DFS at finite temperature for the 
Markovian case [see Fig. U^b)]. We have also checked 
numerically that non-Markovian oscillations of entangle- 
ment gradually diminishes as one increases N defined in 
(p~0|) . for fixed values of other parameters. 



dynamics. We also observe that a non-Markovian deco- 
herence free state (for example, \4>2}) remains decoher- 
ence free in the Markovian regime, but all the Marko- 
vian decoherence free states (for example, |0i)) are not 
necessarily decoherence free in the non-Markovian do- 
main [FigJIJd)]. Finally, we extend our result for the 
finite-temperature case where we see the non-Markovian 
entanglement oscillations gradually decreases as one in- 
creases the temperature. We found that the Markovian 
decoherence- free states remain invariant under finite tem- 
peratures [FigfJJb)]. We also see from FigfJ]that the fi- 
nite temperature of the bath accelerates the phenomenon 
of entanglement sudden death for the non-decoherence- 
free entangled initial states. Interestingly, the state ^2) 
is found to be decoherence free both in the Markov and 
non-Markov regime, and is also robust against finite bath 
temperatures. There were considerable number of papers 
in recent literature |28|, [29|, [3l| dealing with the squeezed 
reservoir. Also, there were proposals J32H361 l38l I39J for 
physical realizations of environments that mimic or gen- 
erate Markovian squeezed bath. One can suggest in fu- 
ture experimental schemes to test and observe various 
important issues related to non-Markovian dynamics by 
physically engineering the environment. Hence, it will be 
interesting to see if our prediction of the non-Markovian 
effects can be verified in actual experiments. 



V. CONCLUSION 

In summary, we study the non-Markovian entangle- 
ment dynamics of two qubits in a common squeezed 
bath. We consider the initial two-qubit states which 
are very close to (as well as far from) the Markovian 
DFS for this system. We see (FigfTJ and FigJ5]) multi- 
ple cycles of entanglement sudden death and revival in 
the non-Markovian case, showing striking difference be- 
tween the Markovian and non-Markovian entanglement 
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